Calculation of some properties of the vacuum 
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In this article, we calculate the dressed quark propagator with the flat bottom potential 
in the framework of the rain-bow Schwinger-Dyson equation, which is determined by mean field 
approximation of the global colour model lagrangian. The dressed quark propagator exhibits a 
dynamical symmetry breaking phenomenon and gives a constituent quark mass about 392 MeV, 
which is close to the value of commonly used constituent quark mass in the chiral quark model. 
Then based on the dressed quark propagator, we calculate some properties of the vacuum, such as 
quark condensate, mixed quark condensate g s {0\qG q\Q) , four quark condensate {0\qTqqVq\0) , 
tensor, ir vacuum susceptibilities. The numerical results are compatible with the values of other 
theoretical approaches. 
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I. INTRODUCTION 



To investigate the long distance properties of quantum chromodynamics (QCD), many approaches are proposed, 
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C ' i such as lattice gauge theory, QCD sum rules, chiral perturbation theory and phcnomcnological quark models. Each 
I i of these approaches has both advantages and disadvantages. For example, lattice calculations are rigorous from 
the point of view of QCD, but they suffer from lattice artifacts and uncertainties connected with the necessary 
, extrapolation to the physical quark masses. The QCD sum rules approach, introduced by Shifman, Vainshtein and 
Zakharov [1,2] in 1979, tries to incorporate nonperturbative elements of full QCD. As the starting point, operator 
product expansion (OPE) method was used to expand the time ordered currents into a series of quark and gluon 
condensates which parameterize the long distance properties, while the short distance effects arc incorporated in the 

Oh, 

QJ , Wilson coefficients. However, these elements (condensates) can not yet be calculated rigorously in QCD, and they are 
(— I \ 

determined from experimental data in the one or other way (For example, the quark condensate can be determined 

! 

. £h [ from experimental data with the help of Gell-Mann-Oakes-Renner relation [1,3].); moreover, although for medium and 
asymptotic momentum transfers the OPE method can be applied for form factors and the moments of wave functions 
[4-6], at low momentum transfer, the standard OPE method cannot be consistently applied as pointed out in the 
early work on photon couplings at low momentum for the nucleon magnetic moments [7,8]. In Ref. [7], the problem was 
solved by using a two point correlator in an external electromagnetic field, with vacuum susceptibilities introduced 
as parameters for nonperturbative propagation in the external field. In fact, the existence and magnitude of the 
vacuum susceptibilities are by themselves an important physical information about the structure of the vacuum and 
have many applications [9-12]. As nonperturbative vacuum properties, the susceptibilities can be introduced for both 
small and large momentum transfers in the external fields . In Ref. [8] , with the long distance effects treated by bilocal 
power corrections and assumption of p meson dominance, the authors circumvent the problem using a three point 
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f The OPE (Wilson expansion) is valid as long as there is a clear distinction between short distances (oc i) which determine 
the Wilson coefficients and long distances (oc -\, it ~ Aqcd) which govern the condensates. When the transformed momentum 
Q 2 is small, the standard expansion in -J^ would lead to a divergent result. 
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formulation. For detailed discussion of the relationship between three point and two point external field treatments 
and the origin of the susceptibilities, one can see Ref. [13]. The two point method, however, has two main problems: it 
cannot be used to extend the coupling to medium and high momentum transfer and there are additional parameters to 
be determined, the vacuum susceptibilities. In Ref. [14], by comparing terms appearing in the two point external field 
expression with those in hybrid expansion of the three point function, the authors obtain a relationship between the 
nonperturbative elements in the two methods. From these relationships, one can express the induced susceptibilities of 
the two point method in terms of well-defined four quark vacuum matrix elements. These susceptibilities may play an 
important role in the QCD sum rules approach. In particular, the strong and parity-violating pion-nucleon coupling 
depends crucially on the tt vacuum susceptibility (x"a), while the tensor vacuum susceptibilities are relevant for the 
determination of the tensor charge of the nucleon, which is connected through deep-inelastic sum rules to the leading- 
twist nucleon structure functions and transversity distribution [15-17]. Although in principle the pion susceptibility 
can be estimated using PCAC, in practice there are inconsistencies with the values needed phenomenologically [14,12]. 
On the other hand, as it has been pointed out in Ref. [18], that treatment of the vacuum tensor susceptibility is subtle 
and different treatments can lead to very different results for the tensor charge. Therefore, it is interesting to calculate 
these vacuum condensates and vacuum susceptibilities within the framework of coupled flat-bottom potential (FBP) 
model [19] and global colour model (GCM) [20-23]. 

The global colour model (GCM), a NJL-like theory [24], has provided a lot of successful descriptions of the long 
distance properties of strong interaction and the QCD vacuum as well as hadronic phenomena at low energy based 
on the theoretical foundation that the quark propagator contains valuable information about the nonperturbative 
properties of QCD [20-23]. The flat-bottom potential is a sum of Yukawa potentials, which not only satisfies gauge 
invariance, chiral invariance and fully relativistic covariance, but also suppresses the singular point which the Yukawa 
potential has. It works well in understanding the meson structure, such as electromagnetic form factor, radius, decay 
constant [19]. 

In this article, we combine the GCM with the FBP to calculate the quark condensate, g s {Q\qG q\Q) , 
(0\qTqqTq\0) , tensor, pion vacuum susceptibilities in the framework of the rain-bow Schwinger-Dyson (SD) equa- 
tion. The article is arranged as follows: in section 2, we brief out line the GCM and in section 3, introduce the FBP, 
SD equation; in section 4, we calculate the vacuum condensates and vacuum susceptibilities; in section 5, conclusion 
and discussion. 



Here we brief out line the main skeleton of the GCM. The global colour model, based upon an effective quark- 
quark interaction which is approximated by the effective flat bottom potential in this article, is defined through a 
truncation of QCD lagrangian while maintaining all global symmetries. 

The generating function for QCD in Euclidean space can be written as 



II. GLOBAL COLOUR MODEL 




(1) 



where 




and G^ v = d^A^ — d v A^ + gf abc ' A h ^A c v . We leave the gauge fixing term, the ghost field term and its integration 
measure to be understood. Here we introduce 
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e W[J] = J VA J ' **{-\G%G%+JlAl) _ (2) 

So the generating function can be written as 

Z[fj,rj\ = J VqVqe- f d4xm ~<- d+m)g - rlq - m] e W[J] . (3) 

The functional W[J] can be formally expanded in the current J": 

W[J] = \J d i xd i yJ«(x)D a Jl(x,y)4(y) + ± J J^J b v J c p D^ p + ■ ■ ■ . (4) 

The GCM is defined through a truncation of the functional W[J] in which the higher order n(> 3)-point functions 
are neglected, and only the gluon 2-point function D^(x,y) is retained. This model maintains all global symmetries 
of QCD and permits a 1/N C expansion, however local SU(3) gauge invariance is lost by the truncation. 
Using the functional integration approach, the generating function of the truncation is given by 

Z% vhcu = J VqVqcxp ( - J rf 4 xg( 7 • 8 + m)q - 9 - J d^yj^D^x - y)j b v (y) + J d 4 x(fjq + q V )), (5) 

where j^(x) — q{x)j^^-q(x) is the quark color current. For simplicity we use a Feynman-like gauge D ab v {x — y) = 
Sfj,iyS ab D(x — y). Performing the standard bosonization procedure, (For detailed discussion of this topic, see Ref. [23]) 
we obtain the resulting expression for the generating function in terms of the bilocal field integration 

Z GC m = J VBe~ s ^\ (6) 

where the action is given by 

S[B] = -TrLnIG- 1 ] + f d 4 xd 4 y ^^^ + J d 4 xd 4 yfj(x)G(x, y; [B e }) V (y), (7) 

and the inverse quark's Green function G _1 is defined as 

G-\x,y) = 1 -dS(x-y)+A e B e (x,y). (8) 
Here the quantity A e arises from Fierz reordering of the current-current interaction term in Eq.(5), 

^-in^-mj — (lti~^)ij{ln~^)mn (9) 

and is the direct product of Dirac, flavor SU(3) and color matrices, 

A e = i(lz^75, ^=7m> -jf ^75) ® (^=1f, ^=A?0 ® (^l c , -J=A«). (10) 



The vacuum configurations are defined by minimizing the bilocal action: SS g ^ 
field approximation, we obtain 



= 0, following this mean 

fj.ri— 



B 6 (x - y) = 9 2 D(x - y)tr[K 9 G {x - y)]. (11) 

These configurations provide self-energy dressing of the quarks through the definition S(p) = k e B®{p) = ij ■ p[A(p 2 ) — 
1] + B(p 2 ). In fact, equation (11) is the rain-bow SD equation. In terms of A and B, the quark's Green function at 
Bq is given by 

Gn(x y) - Gn(x - v) - f ^ ■ P A (P 2 ) + B(p 2 ) (x _ y) 

G (x lV )-G (x y)-J ^ p2A2{p2) + B2{p2) e . (U) 
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With the GCM generating function Eq.(6), it is straightforward to calculate the vacuum expectation value of any 
quark operator of the form 

O n ee (fcAW^JfeAS.fe) • • • (^Ag^-J, (13) 

in the mean field vacuum. Here the AW represents an operator in Dirac, flavor and color space. 

Taking the appropriate number of derivatives with respect to the external sources r\i and ffj of Eq. (6) and putting 
i]i = fjj = [25], we obtain 

(0| : O n : |0) = (-)" ^[(-) p A« • • • a£> b (G Q ) hip(1) ■ ■ ■ (G ) jnipi J, (14) 
p 

where p stands for a permutation of the n indices. From Eq.(14), we can easily obtain the expression for the nonlocal 
four quark condensate, which is another important vacuum condensate in the QCD sum rules approach beside quark 
condensate. 

(0| : q(x)A^q(x)q(y)A^q(y) : |0) - -tr jC [G (y, x)A^G (x, y)A^} + tr^ c [G (x, x)A^}tr, c [G (y, y)A {2) }. (15) 



III. FLAT BOTTOM POTENTIAL AND SCHWINGER-DYSON EQUATION 

The phenomenological FBP is a sum of Yukawa potentials which is an analogous to the exchange of a series of 
particles and ghosts with different masses, 

n 

where N stands for the minimum value of the mass, p is their mass difference, and aj is their relative coupling 
constant. In its three dimensional form, the FBP takes the following form: 

™ e -(N+jp)r 

V(r) = -J2a j . (17) 

3=0 



In order to suppress the singular point at r = 0, we take the following conditions: 

V(0) = constant, 

dv(o) = d 2 v(o) = = rf-y(o) = Q 

dr dr 2 dr n 

So we can determine aj by solve the following equations, inferred from the flat bottom condition Eq.(18), 

n n n n 

5>i=0, ^a,(iV + j» = y(0), ^a,(7V + . ?P ) 2 = 0, ••• aj (N + jp) n = 0. (19) 

3=0 ]=0 j=0 j=0 

As in previous literature [19], n is set to be 9. In Rcf. [19] (Phys. Rev. D47 (1993) 2098), the authors solve the 
Bethe-Salpetcr equation based on expansion in terms of Gegenbauer polynomials with FBP, then use the Bethe- 
Salpctcr amplitudes as input to calculate the 7r meson electromagnetic form factor. Among many other values, the 
phenomenological value 9 can give stable solution and best fit to experimental data* . 
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In the rainbow approximation, the SD equation takes the following form: 




(20) 



where 



S-'ip) = A(p 2 ) 7 • p - B(p 2 ) ee Aip 2 )^ ■ p - m(p 2 )} 
= -g^G(k 2 ), 



(21) 
(22) 



and m stands for an explicit quark mass-breaking term. With the explicit small mass term, we can preclude the 
zero solution for B(p) and in fact there indeed exist a bare current quark mass. In order to keep consistently with 
section two, here we take Feynman gauge. This dressing comprises the notation of constituent quarks by providing a 
mass m(p 2 ) = B(p 2 )/A(p 2 ), which is corresponding to dynamical symmetry breaking. Because the form of the gluon 
propagator g 2 G(p) in the infrared region is unknown, one often uses model forms as input in the previous studies of 
the rainbow SD equation [20,22,23,36]. 

As a nonperturbative phenomenological potential, the FBP may have intrinsic connections with other nonpertur- 
bative phenomena, such as confinement, instanton, vacuum condensate and renormalon, glueball. In this article, we 
concern mainly the quark condensates, while the vacuum condensates play an important role in determining the prop- 
erties of the mesons, glueballs and their mixing as indicated by the low energy theorem [28]. However, the realization 
of gluonic freedom in hadronic models is more ambiguous and difficulty then that for quarks, the glueball spectra 
obtained from two main nonperturbative methods (lattice QCD and QCD sum rules ) have some contradictions [30]. 
Here we just replace the gluon propagator with FBP and suppose the nonperturbative effects are incorporated in it, 
the possible effects on glueball will be our next work. § Although there may be debates against doing that, the quark 
propagator obtained with the FBP does lead to satisfactory dynamical symmetry breaking phenomenon and quark 
confinement based on the theoretical foundation that there no singularities on the real timelike p 2 axial [29] . 

In this article, we assume that a Wick rotation to Euclidean variables is allowed, and perform a rotation analyti- 
cally continuing p and k into the Euclidean region where them can be denoted by p and k, respectively. The Euclidean 
SD equation can be projected into two coupled integral equations for A(p 2 ) and B(p 2 ). For simplicity, here ignore 
the bar on p and k. Numerical values for A(p 2 ), B(p 2 ) and m(p 2 ) are shown in Fig.[l] . 



In this section, we calculate the quark condensate (0|qq|0), the mixed quark gluon condensate g s (0\q~G q\0) , 
the four quark condensate (0\qTqqTq\0} and tensor, pion vacuum susceptibilities. The numerical results are compared 
with other theoretical works. 

Here we take a short digression to discussing gauge invariance. The propagation of a colored quark in the 
background gluon field results in a phase factor given by the exponential of a line integral, 



''For example, we can obtain the gluon condensate from the nonperturbative gluon propagator (the effective potential), then 
based on the gluon condensate and low energy theorem, we can study the lowest lying glueball mass, their mixing, and so on. 
There is another way to study the flat bottom potential to glueball spectrum through the Bethe-Salpeter equation. 



IV. VACUUM CONDENSATES AND VACUUM SUSCEPTIBILITIES 




(23) 
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where P is a path order operator, stands for that the color matrices should be ordered along the path C connecting 
the point x and y. The gauge dependence of the quark propagator can be factored out as an " A-depcndcnt" phase 
before taking the vacuum expectation value. In operator expanding of the correlators of colorless currents, these 
phase factors can be canceled out with each other, one can facilitate the calculation by imposing Fock-Schwingcr 
gauge condition [31,32], 

{x»-x»)A ll {x) = Q, (24) 

or in the other word, let P = 1. To gauge dependent quantities, (for example, (0\q(x)q(0)\0)) , we can insert phase 
factors P to produce the corresponding gauge invariant ones ( (0\q(x)P(x, 0; A, C)q(0)\0)). The gauge invariant 
quantity (0| : q(x)P(x, 0; A, C)q(0) : |0) calculated in a special gauge, for example, Fock-Schwinger gauge, will not 
change its value. In the following, the gauge dependent quantities can be taken as their corresponding invariant 
ones in the special gauge, Fock-Schwingcr gauge. Furthermore, although the gauge dependent quantities, such as 
(0| : q(x)q(0) : |0), are dependent on the special gauge we choose, they are important imformations in describing the 
nonperturbative QCD vacuum. In a sense, the gauge dependence is fictitious, since the vacuum field, after averaging, 
is clearly invariant under transformations. 



A. Vacuum Condensates 



The quark propagator is defined as 

S(x) = (0\T[q(x)q(0)}\0)- (25) 

where q(x) is the quark field and T the time-ordering operator. For the physical vacuum consists of both perturbative 
and nonperturbative parts, so the quark propagator S(x) can be divided into a perturbative and a nonperturbative 
part as the following: 

S(x) = S PT (x) + S NP (x). (26) 

In the nonperturbative vacuum, the normal-ordered product Snp(x) does not vanish. For short distance, the OPE 
for the scalar part of Snp(x) gives 

<0| : q(x)q(0) : |0> = <0| : g(0)«(0) : |0) - ^«0| : g(0)a • G(0)«(0) : |0> + ■ ■ ■ (27) 

in which the local operators of the expansion are the quark condensate, the mixed condensate, and so forth. In Ref. 
[34] the nonlocal condensate is put in the following form: 

(0| : q(x)q(0) : |0) = 5 (a; 2 )(0| : q(0)q(0) : |0), (28) 

where g(x 2 ) is the vacuum non-locality of the nonlocal quark condensate. The nonlocal quark condensate (0| : 
q(x)q(0) : |0) is given then by the scalar part of the Fourier transformed inverse quark propagator, which can be 
inferred from Eq.(14). 

(0\:q(x)q(0):\0) = (-)tr 7 c[G o (x,0)] 



(-4Ay / 

Jo 



" d 4 p B( P 2 ) 



(-) 



o (2tt) 4 p 2 A 2 (p 2 ) + B 2 (p 2 ) 
12 f B(s) J^sx 2 )^ 



16tt 2 



r dss m [2 w) , (29) 

J dSS sA 2 ( s ) + B 2 (s) [2 ^2 J' ^ 
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At x=0 the expression for the local condensate (0| : q(0)q(0) : |0) is recovered, 

<0| : g(0) g (0) : |0> = - ^ ^ <toa _^_ . (30) 

The non-locality g(x 2 ) can be obtained immediately by dividing Eq.(29) through Eq.(30). Our analysis ignores 
effects from hard gluonic radiative correction to the condensates which are connected to a possible change of the 
renormalization scale /i at which the condensates are defined. As in Ref. [36], /i is taken to be 1 GeV 2 . Final we 
obtain the nonlocal quark condensate which shown in Fig. [2] and local quark condensate which shown in Table 1. 
From Fig. [2], we can see that the nonlocality is compatible with other theoretical calculations, including instanton 
liquid model [37], effective quark-quark interaction model [34], GCM [38] , dipole fit approach [14], and other potential 
model [41]. The curves obtained in Ref. [34] and in Ref. [14] are similar to each other, here we choose the curve in 
Ref. [34] to compare with our result. For simplicity, in Fig. 2, we list three curves. 

In [44] , we calculate the nonlocality of the quark condensate; in calculation, wc find that there are other parameters 
can give correct nonlocality curve. In the present work, we re-fit the parameters to incorporate the mixed condensate 
and the vacuum susceptibilities. Furthermore, in the present work, we change the coefficient of FBP inserted in the 
SD equation by a factor n, the present formulation is better than the old one [44,19]. 

By means of Eq.(15), one can calculate all kinds of nonlocal four quark condensates . If wc takeA^ =A< 2 )= 7m ^ 

then 

(0\:q(xh^q(x)q(0)l^q(0)--\0) 

= (_) f f fP 4 d4 l c i*-(v-«) La B(p 2 ) B(q 2 ) 
Jo Jo 



/ J (2^) 4 (2tt) 



A 2 (p 2 )p 2 + B 2 (p 2 ) A 2 (q 2 )q 2 + B 2 (q 2 ) 



2 A(p 2 ) A(q 2 ) 

+ X A 2 (p 2 )p 2 + B 2 (p 2 )A 2 {q 2 )q 2 + B 2 (q 2 ) P ' q 



(31) 



Similarly, at x=0 wc obtain the expression for the local four quark condensate (0| : q^n^-qq^^^-q ■ |0), 

(01 : 57,f m4 , : 10) - ^ ^ffi'g^ l' = H^h ft : VP- <*>) 

As far as the mixed condensate g s (qG^a^ v q) is concerned , the evaluation is somewhat lengthy, one can use the 
method described by Ref. [36] to obtain the mixed condensate in Euclidean space, 

In Table 1. we display our results for (qq) and g s {q~G ' q) and compare them with the corresponding values 
obtained from other theoretical approaches, such as QCD sum rules [26], quenched lattice QCD [33], the instanton 
liquid model [35] , the model of confining gluon propagator for GCM [36] and other potential model [39-41] . From 
Table 1, we can see that the values are compatible with other theoretical works. 

B. Vacuum Susceptibilities 

In the external field of QCD sum rule two-point method, one often encounters the quark propagator in the 
presence of a external current J r (y) — q(y)Tq(y)(T stands for the appropriate combination of Dirac , flavor and 
colour matrices). 

S c a f(x) = {0\T[&{x)4{oW)jr = S c :^ PT {x) + S c a ^ NP (x), (34) 
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where r ' PT (x) is the quark propagator coupled perturbatively to the current and S^ c r ' NP (x) is the nonperturba- 
tive quark propagator in the presence of the external current J T (To be explicitly, here we keep all the indexes). 
The vacuum susceptibility x r m the QCD sum rule two-point external field treatment can be defined as [f 4] 

S C a C '^ NP {x) = <0| : q c a (x)4(0) : |0) 7 r = -±T a p6 cc , X T H(x)(0\ : g(0)g(0) : |0>, (35) 

where the phenomenological function H(x) represents the nonlocality of the two quark nonlocal condensate. Obviously 
H(0)=1. 

By comparing terms appearing in the two point external field expression with those in hybrid expansion of the 
three point function, one can obtain a relationship between the nonperturbative elements in the two methods. From 
these relationship, one can express the induced susceptibilities of the two point method in terms of well-defined four 
quark vacuum matrix element [14]. 

The presence of external field implies that 5*^ r (x) is evaluated with an additional term AL = —J r -<pr added 
to the usual QCD Lagrangian, where (f>r is the value of external field. In three point method of QCD sum rule, if one 
takes only a linear external field approximation, the S^ g T ' NP (x) in Euclidean space is given by 



~,cc'T,NP ( 
5 q/3 



\x) = J d 4 y e-^<0| : q c a {x)q{y)Tq{y)4 (0) : |0>. (36) 

Using Eq.(35) and Eq.(36) we obtain 

-^S cc ,T af}X r (0\ : 9(0) g (0) : |0> = J d 4 y e"*™ <0| : <f a (0)q{y)r q {y)$ (0) : |0). (37) 

Multiplying Eq.(37) by Tp a S cc >, we get the result: 

X F a = -~^y / d 4 y e-^(0\ : q(0)rq(0)q(y)r q (y) : |0). (38) 

In the case of tensor current (r = u^), according to Eq.(14), we calculate the tensor matrix element explicitly, 

(0\q{0)(T lu/ q{0)q(y)a IM ,q{y)\0) 

f d 4 p -i-y ■ pA(p 2 ) + B(p 2 ) f d 4 fc _ -i-y • kA{k 2 ) + B{k 2 ) 



// 



(2tt) 4 ^ A(p 2 )p 2 + B{p 2 ) 7L J (2tt) 4 ^ A(k 2 )k 2 + B(k 2 ) 
d 4 p d 4 fc __ i(j ,- k ).„^ L • PA{ P 2 ) + B(p 2 ) _ -ij-kA(k 2 , - B(^) 



(2tt) 4 (2tt) 4 



L 7C 



'"" A{p 2 )p 2 + B(p 2 ) ^ A(k 2 )k 2 + B(k 2 ) 



A8N I I *E-*Le-*0-»-y ^ (39) 

c/ ' (2tt) 4 (2^) 4 A(p 2 )p 2 + B(p 2 ) A(k 2 )k 2 + B{k 2 Y K ' 



where N c = 3 is the number of colors. 
In the limit q — > 0, we obtain 



r 

X a — 3 / sds 
Jo 



B(s) 



2 



= 0.0757GeV . (40) 



A 2 s + B 2 (s)_ 

From Eq.(36) we can see that the vacuum susceptibility originates from the nonlocal four quark condensate 
contribution. In Ref. [43], the authors get the same conclusion based on completely different viewpoint of duality. 

The vacuum susceptibility x T a defined in Eq.(38) (same as Ref. [39]) has an opposite sign and a factor Air 2 larger 
than the definition in Refs. [15,18,43]. In order to compare with those estimations, we make a redefinition 

X T a -> = -0.0019GeU 2 . (41) 
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In Table 2, we compare the result with the values of other theoretical estimations of the vacuum tensor suscep- 
tibility. 

In the following, we calculate the 7r vacuum susceptibility which is crucial to determine the strong and party- 
violating pion-nuclcon coupling. In the case of pseudoscalar current, we take AW = A< 2 ) = 7 5 in Eq.(15), and get the 
7r vacuum susceptibility x" a - 



X^a = 12-k 2 / d 4 y tr 7 
1 



\dStr, [J ^ 



d A p -i-y ■ pA + B(p 2 ) 



Jo 



sds- 



) 4 A 2 p 2 + B 2 (p 2 ) 
= 1.06 GeV 2 , 



15 J W 



A -i>y. q A + B(q 2 ) lq . y 
) 4 A 2 q 2 + B 2 (q 2 ) 



75 



(42) 



A 2 s + B 2 {s) 

which is below the range x* a — (1-7 — 3.0) GeV 2 obtained within a phenomenological approach [14]. In Rcf. [14], 
the author made a simple estimation of the value based on the space time structure of the nonlocal quark condensate 
extracted from experimental data on sea-quark distributions but with a crude assumption of vacuum saturation for 
the intermediate states. 



V. CONCLUSION AND DISCUSSION 



In this article, we calculate the dressed quark propagator with the FBP in the framework of the rain-bow SD 
equation, which is determined by mean field approximation of the GCM lagrangian. The dressed quark propagator ex- 
hibits a dynamical symmetry breaking phenomenon and gives a constituent quark mass about 392 MeV, which is close 
to the value of commonly used constituent quark mass 350 MeV in the chiral quark model. Then based on the dressed 
quark propagator, we obtain the quark condensate (0|qq|0), the mixed quark gluon condensate g s (0\q~G q\0) , the 
four quark condensate (0\qTqqTq\0) and tensor, pion vacuum susceptibilities at the mean field level, which are com- 
pared with other theoretical calculations. The values of the nonlocality of quark condensate, local quark condensate, 
mixed quark condensate are compatible with other theoretical results. While the value of tt susceptibility is below 
that estimated in Rcf. [14], which is based on the space-time structure of the nonlocal quark condensate and a crude 
assumption of vacuum saturation for the intermediate states. In Table 3, we also compare with other theoretical 
results. In Table 2, we can see that the numerical value of the tensor vacuum susceptibility varies with theoretical 
approach. It has been pointed out in Ref. [15] that experimental measurement of the tensor charge of the nucleon is 
possible, and one might be able to test the theoretical predictions of the tensor susceptibility in the future. 

In solving the rain-bow SD equation, the parameters are taken as Aqcd — 200MeV,N = 3.0Aqcd, V(0) = 
— 12.0Aqcd, p = 2.0Aqc_d, m u = ma = 8MeV and the large momentum cut-off L = 630(Aqcd)- 

In calculation, the coupled integral equations for the quark propagator functions A(p 2 ) and B(p 2 ) are solved 
numerically by simultaneous iterations. The iterations converge rapidly to a unique stable solution of propagator 
functions and independent the initial guesses. The propagator functions A(p 2 ) and B(p 2 ) are shown in Fig.[l] , at 
small p 2 , A(p 2 ) differs from the value 1 appreciably, while it tends to 1 for large p 2 . We find that at small p 2 , m(p 2 ) 
is greatly re-normalized, while at large p 2 , it takes asymptotic behaviour. For u and d quark, m(0) = 392MeV, which 
is close to the constituent quark masses, the connection of m(j>) to constituent masses is somewhat less direct and is 
precise only for heavy quarks. For heavy quarks, m constituent (p) = m(p = 2m constituent (p)) , for light quarks , it only 
makes a crude estimation [42]. At about p = lGeV, the mass function grows rapidly as the momentum decreases, 
that is an indication of dynamical symmetry breaking. 

The results show that the phenomenological FBP works well with the rainbow SD equation in calculating nonlocal 
quark vacuum condensate, mixed condensate and vacuum susceptibilities. This can be extend to other system, such 
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as qq system, qQ system, vector meson and qqq system or by applying the results to the calculation of quantities and 
processes requiring detailed knowledge of the quark propagator . 
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Fig. 1 SD functions 
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Table 1: The values of {0\qq\0), g(0\qaGq\0). 



Reference 


-(qq)* 


-g(qaGq)i 




[MeV] 


[MeV] 


[26] 


210-230 


375-395 


[33] 


225 


402-429 


[35] 


272 


490 


[36] 


150-180 


400-460 


[39] 


217 


429 


[40] 


178 


456 


[41] 


171 


447 


This Work 


182 


422 



1 



Table 2: The values of x T(l 



Reference 


T 

X a 




[GeV] 


[15] 


0.002 


[18] 


-0.008 


[27] 


0.009-0.017 


[43] 


-0.0055 


[40] 


-0.0014 


[39] 


-0.0058 


[41] 


-0.0018 


This Work 


-0.0019 



Table 3: The values of x^ a 



Reference 






[GeV] 


[14] 


1.7-3.0 


[11] 


1.88 


[39] 


2.56 


[41] 


0.68 


This work 


1.06 



2 



